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We review finite-time future singularities in modified gravity. We reconstruct an explicit model 
of modified gravity realizing a crossing of the phantom divide and show that the Big Rip singular- 
ity appears in the modified gravitational theory. It is also demonstrated that the (finite-time) Big 
Rip singularity in the modified gravity is transformed into the infinite-time singularity in the corre- 
sponding scalar field theory obtained through the conformal transformation. Furthermore, we study 
several models of modified gravity which produce accelerating cosmologies ending at the finite-time 
future singularities of all four known types. 



I. INTRODUCTION 



X 



On 
O 

, 

O^l , Recent observations confirmed that the current expansion of the universe is accelerating. There are two approaches 
to explain the current accelerated expansion of the universe. One is to introduce some unknown matter, which is 
or called "dark energy" in the framework of general relativity. The other is to modify the gravitational theory, e.g., to 

■ study the action described by an arbitrary function of the scalar curvature R. This is called ll F(R) gravity" , where 
F(R) is an arbitrary function of the scalar curvature R (for reviews, sec M). 

According to the recent various observational data, there exists the possibility that the effective equation of state 
(EoS) parameter, which is the ratio of the effective pressure of the universe to the effective energy density of it, evolves 

1 ^ ', from larger than —1 (non-phantom phase) to less than —1 (phantom one, in which superacceleration is realized), 
^3 ' namely, crosses —1 (the phantom divide) currently or in near future. Various attempts to realize the crossing of the 

JL , phantom divide have been made in the framework of general relativity. Recently, a crossing of the phantom divide 
'■ in modified gravity has also been investigated [l], H, 0,0- Moreover, it is known that modified gravity may lead to 
\ the effective phantom/quintessence phase while the phantom/quintessence-dominated universe may end up with 
l— ~~ '■ finite-time future singularities, which can be categorized into four types 0]. 

In thepresent article, we review finite-time future singularities in modified gravity. Following the considerations 
in Ref . [5| , we reconstruct an explicit model of modified gravity realizing a crossing of the phantom divide by using 
I/*-) ' the reconstruction method proposed in Refs. 0, [1| • We show that the Big Rip singularity appears in this modified 
, gravitational theory, whereas that the (finite-time) Big Rip singularity in the modified gravity is transformed into 
• the infinite-time singularity in the corresponding scalar field theory. Next, following the investigations in Ref. Q, we 
explore several examples of F(R) gravity which predict the accelerating cosmological solutions ending at the finite- 

■ time future singularities. It is demonstrated that not only the Big Rip but other three types of the finite-time future 
' singularities may appear. 

0^ . This article is organized as follows. In Sec. II we explain the reconstruction method of modified gravity @, [1]. 
Using this method, we reconstruct an explicit model of modified gravity in which a crossing of the phantom divide 
, can be realized. We also consider the corresponding scalar field theory. In Sec. Ill we present several models of F(R) 
gravity which predict accelerating cosmologies ending at the finite-time future singularities by using the reconstruction 
method. Finally, summary is given in Sec. IV. We use units in which ks = c = H = 1 and denote the gravitational 
constant 8ttG by k 2 , so that n 2 = 8tt/M p1 2 , where Af P1 = G~ 1/2 = 1.2 x 10 19 GcV is the Planck mass. 



II. MODIFIED GRAVITATIONAL THEORY REALIZING A CROSSING OF THE PHANTOM DIVIDE 



We investigate modified gravity in which a crossing of the phantom divide can be realized by using the reconstruction 
method. 



A. Reconstruction of modified gravity 



First, we review the reconstruction method of modified gravity proposed in Refs. @, 1M| (for the related study of 
reconstruction in F(R) gravity, see (Toj). 
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The action of F(R) gravity with general matter is given by 

'F(R) 



I 



S = Id X\J — g 



2k 2 



(2.1) 



where g is the determinant of the metric tensor g^ v and £ ma tter is the matter Lagrangian. 

The action (|2.ip can be rewritten to the following form by using proper functions P(<p) an d Q(4>) of a scalar field 



S = J d 4 .T^|^2 [P(4>)R + Q{4>)\ + Anatterj • 



(2.2) 



The scalar field <p may be regarded as an auxiliary scalar field because <p has no kinetic term. Taking the variation of 
the action (|2.2p with respect to <j), we obtain 

dip dip 

which may be solved with respect to <p as <p = <p(R). Substituting <p = <p(R) into the action (|2.2I) . we find that the 
expression of F(R) in the action of F(R) gravity in Eq. (|2.ip is given by 

F(R)=P(<f>(R))R + Q(c/>(R)). (2.4) 

Taking the variation of the action (|2.2p with respect to the metric g^ v , we find that the field equation of modified 
gravity is given by 

hg^ [PWR + Q{4>)] - R^P{4>) - 9^P(<P) + V M V„P(0) + k 2 T$^ = , (2.5) 

where V M is the covariant derivative operator associated with g^, □ = g^V ^ v is the covariant d'Alembertian for 
a scalar field, and j^™ atter ) j s t ne contribution to the matter energy-momentum tensor. 

We assume the flat Friedmann- Robertson- Walker (FRW) space-time with the metric ds = —dt + a 2 {t)dx 2 , where 
a(t) is the scale factor. In this background, the (fJ>,v) = (0,0) component and the trace part of the (fJ>,v) = 
component of Eq. ()2.5|) . where i and j run from 1 to 3, read 

-&H 2 P(cp{t)) - WW) - 6H dP( f (f)) + 2k 2 p = , (2.6) 

dt 

and 2 d 2 PW)) +iH dP(m) + U il + &H 2\ Pm ))+Q(cb(t)) + 2 K 2 p = 0, (2.7) 

dt A dt \ J 

respectively, where H = a/a is the Hubble parameter. Here, p and p are the sum of the energy density and pressure 
of matters with a constant EoS parameter Wi, respectively, where i denotes some component of the matters. 
Eliminating Q(<p) from Eqs. (|2.6p and (|2.7p . we obtain 

d 2 P(6(t)) dPU(t)) ■ , , s , 9 , , s 

~ ~ H dt + 2HP ^ + (P + P) = ■ (2-8) 

We note that the scalar field </> may be taken as <p = t because <p can be redefined properly. 

We now consider that a(t) is described as a(t) = acxp (g(t)), where a is a constant and g(t) is a proper function. 
In this case, Eq. (|2.8p is reduced to 



d 2 P{4>) d~g{4>) dP{(j)) ^ d 2 g((p) jo,. , , an+to,) to/-,, n to a^ 

~1^ 2 cup — If ~df ^ w w 5 = ' 

where is a constant and we have used H = dg(<p)/ (dip). Moreover, it follows from Eq. (|2.6p that Q{<p) is given by 



W) 



d~g{<P) 



P{4>) + 2k 2 Pi"- 3(1+m) ^P [-3 (1 + m) m • (2-10) 



Hence, if we obtain the solution of Eq. (|2.9p with respect to P(<p), then we can find Q(<p). Consequentl y, u sing 
Eq. (|2.4p . we can reconstruct F(R) gravity for any cosmology expressed by a(t) = aexp (g(t)). In Refs. 111 El, 
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specific models unifying the sequence: the early-time acceleration, radiation/mattcr-dominatcd stage and dark energy 
epoch have been constructed. 

Next, using the above method, we reconstruct an explicit model in which a crossing of the phantom divide can be 
realized. We start with Eq. (|2.9p without matter: 



d 2 P(<f>) dg(^)dP(cb) , d 2 g(^) 



P(0) = . 



(2.11) 



By redefining P(<j>) as P(4>) 



p{4>); Eq. (|2.11[) is rewritten to 



1 d 2 p{cj>) 
p(4>) dip 2 



25e 9(«/10 



d 2 (e 



-g(4>)/w\ 



^7+l 



We explore the following model: g{<j>) = —10 In (4>/to) 7 — C (4>/to 
is the present time. In this case, Eq. (j2~T2)l is reduced to (1 /§((/>)) [d 2 p{(j))/ {d<p 2 )} = 257(7 + 1)/0 2 , which can be solved 



(2.12) 



where C and 7 are positive constants, and to 



as p(4>) = p~j r 4> l3+ + p~(pP~ . Here, p± are arbitrary constants and j3± are given by (3± = 1 ± \/l + 1007(7 + 1) /2- 

From the above expression of g(4>), we find that g{4>) diverges at finite </> when <f> = t s = toC^ 1 ^ 2l+1 \ which tells that 
there could be the Big Rip singularity at t = t s . We only need to consider the period < t < t s because g(<f>) should 
be real number. In this case, the Hubble rate H(t) is given by 



H(t) 



dm 



to 



-7-1 



to 



to 



C 



to 



7+1' 



(2.13) 



where it is taken cf> = t. 

In the flat FRW background, even for modified gravity described by the action (|2.1j) . the effective energy-density 
and pressure of the universe are given by p e g = 3H 2 / k 2 and p e g = — (2H + 3H 2 ^j /k 2 , respectively The effective 

EoS parameter w e g — p e s/peS is defined as w c r = — 1 — 2H / (3-ff 2 ) [l|. For the case of H(t) in Eq. (|2. 13|) . w c s is 
expressed as w e s = —1 + U(t), where 



U(t) 



2H 
3H 2 



1 

15 



-7 + 47(7 + 1) ( - 



27+1 



+ (7 + j 



2(27+1) 



7 + (7 + 1) ( - 



27+1 



(2.14) 



For the case of Eq. (|2.13[) . the scalar curvature R = 6 ( H + 2H 2 ) is expressed as 



R = 60 



7 (20 7 -l) + 44 7 (7 + 1) 



27+1 



+ ( 7 +l)(20 7 + 21) 



t x 2(27+1)' 



27+I' 



(2.15) 



In deriving Eqs. and (pjB) . we have used t s = t C" 1 / (27+1 ). 

When t — * 0, i.e., t < t s , H(t) behaves as H(t) ~ lOy/t. In this limit, it follows from w cS = -1 - 2H / (3i? 2 ) that 
the effective EoS parameter is given by w e s = —1 + 1/ (157). This behavior is identical with that in the Einstein 
gravity with matter whose EoS parameter is greater than —1. 

On the other hand, when t — ► t s , we find H(t) ~ 10/ (t s — t). In this case, the scale factor is given by a(t) ~ 
a(t a — t) 10 . When t — > t s , therefore, a — > 00, namely, the Big Rip singularity appears. In this limit, the effective 
EoS parameter is given by w c g = —1 — 1/15 = —16/15. This behavior is identical with the case in which there is a 
phantom matter with its EoS parameter being smaller than —1. Thus, we have obtained an explicit model showing 
a crossing of the phantom divide. 

It follows from w e g = — 1 — 2H / (3-ff 2 ) that the effective EoS parameter w e g becomes —1 when H = 0. Solving 

w c ff = —1 with respect to t by using w c s = —1 + U(t), namely, U(t) = 0, we find that the effective EoS parameter 

1/(27+1) 



crosses the phantom divide at t = t c given by t c 



t. 



-2 7 + V47 2 + 7/ (7 + 1) 



From Eq. (|2.14j) . we see 



that when t < t c , U(t) > because 7 > 0. Moreover, the time derivative of U(t) is given by 



dU{t) 
dt 



1 

15 



27 (7+I) (27 + 1)' 



7 +(7+1) 



27+1' 



-1 -3 



27 



27+1' 



(2.16) 
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Eq. (|2.16[) tells that the relation dU(t)/ (dt) < is always satisfied because we only consider the period < t < t s 
as mentioned above. This means that U(t) decreases monotonously. Thus, the value of U(t) evolves from positive to 
negative. From w e g = —1 + U(t), we see that the value of w e g crosses — 1. Once the universe enters the phantom 
phase, it stays in this phase, namely, the value of w e g remains less than — 1, and finally the Big Rip singularity appears 
because U(t) decreases monotonically. 

As a consequence, P(t) is given by P(t) = j(t/t ) 7 /I - {t/t s f 1+1 } Ej=± • Usin S Ec l s - <EHSJ), we obtain 



Q(t.) = -6H {(t./t. Q y /[l -(*/*.) 



When t -> 0, from H(t) ~ lOf/t, we find t ~ ^607 (2O7 - 1) JR. In this limit, it follows from Eqs. dUD that the 
form of F(R) is given by 



F(R) 



^607(207- 1)^ 



-1/2 



1 - 



f^607(20 7 - TjR- 1 / 2 



27+1 



^E 

j'=± 



/ 5 7 -l-/? 3 - 

V 207 - 1 



p J [60 7 (207-l)] / 



ft/ 2 i^-ft/2 



(2.17) 



On the other hand, when t — ► i s , from -ff(t) ~ 10/ (t s — i), we obtain t ~ t s — 3-\/140/i?. In this limit, it follows 
from Eqs. (|2.4p that the form of F(R) is given by 



F(R) 



(J/toY 



1 - {J /ts) 



27+I 



3=± 



1 - 



\t a + {/3 s - 15) R- 1 ' 2 



J 



(2.18) 



where J = t s — Zy/\4QjR. The above modified gravity may be considered as some approximated form of more realistic, 
viable theory. From Eq. (|2.15[) . we see that in the above limit the scalar curvature diverges, and that the expression 
of F(R) in (|2.18j) can be approximately written as 



F(R) 



1 

3a/T40(2 7 + 1) V*o 



1 5 



VJ = ± 



(2.19) 



B. Corresponding scalar field theory 

In this subsection, motivated by the discussion in Ref. [l3|, we consider the corresponding scalar field theory to 
modified gravity realizing a crossing of the phantom divide, which is obtained by making the conformal transformation 
of the modified gravitational theory. (In Ref. [l4|, the relations between scalar field theories and F(R) gravity have 
been studied.) By introducing two scalar fields C and £, we can rewrite the action (|2.1[) to the following form [1[: 




(2.20) 



The form in Eq. (|2 . 20[) is reduced to the original expression in Eq. (|2.ip by using the equation ( = R, which is derived 
by taking variation of the action (|2.20[) with respect to one auxiliary field £. Taking the variation of the form in 
Eq. (|2.20|) with respect to the other auxiliary field £, we obtain £ = F'((), where the prime denotes differentiation 
with respect to (. Substituting this equation into Eq. (|2.20[) and eliminating £ from Eq. (|2.20p . we find 



S 



— (F'(0R + F(C) - F'(C)C) + £ matter 



(2.21) 



This is the action in the Jordan- frame, in which there exists a non-minimal coupling between F'(£) and the scalar 
curvature R. We make the following conformal transformation of the action (|2.2ip : — > = e a g^ Vl where 
e a = F'(£). Here, a is a scalar field and a hat denotes quantities in the Einstein frame, in which the non- minimal 
coupling between F'(() and the scalar curvature R in the first term on the right-hand side of Eq. (|2.21[) disappears. 
Consequently, the action in the Einstein frame is given by [Tol IT(| 



Sk = 



( 



(2.22) 
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where V(a ) = e- ff C(c) - er'^F «(<r)) = (/F'(Q - F(Q/ (F'(C)) 2 and g is the determinant of g^ '. In deriving 
Eqs. (|2.22[) . we have used e a = F'(Q. In addition, ((a) is obtained by solving e a = F'(() with respect to f as 
C = ((f )- Defining f as (p = y // 3/2a/n, the action (|2.22p is reduced to the following form of the canonical scalar field 
theory: 

<Sst = / d ^ Xy 



R 



(2.23) 



Taking the variation of the action ()2.1() with respect to the metric g^ v , we find that the field equation of modified 
gravity is given by F'(R)R^ - (1/2) g flv F{R) + g^UF'{R) - V M V„F'(,R) = K 2 T^ atter) . When there is no matter, 
using the (//, v) = (0, 0) component and the trace part of the (/x, v) = (i, j) component of the above gravitational field 
equation, in the flat FRW background, we obtain 

2HF'(R) + 6 (-4H 2 H + 4H 2 + 3HH + Iff) F"(R) + 36 UhH + Hj F"'{R) = . (2.24) 

We now investigate the case in which F(R) is given by F(R) = Ci M 2 (R /M 2 ) ™, where ci is a dimension- 
less constant and M denotes a mass scale. The form of F(R) in Eq. (|2.19p corresponds to the above expres- 
sion with n = —7/2. In this case, the scale factor a(t) and the scale curvature R are given by [l3| a(t) = 
a (t s - t)(«+ 1 )(2™+ 1 )/(«+ 2 ) an d R = 6n ( n + !)( 2n + i)(4 n + 5)( n + 2)" 2 (t s - t)~ 2 , respectively. It follow from 
dt = ±e a > 2 dt that the relation between the cosmic time in the Einstein frame i and one in the Jordan frame is 
given by t = ^yf=nci{n + 2)" [6n(n + l)(2n + l)(4njJ5)]~ (n+1)/2 M Tl+1 (i s - t) ,l+2 . If n < -2, the limit of t -> t s 
corresponds to that of < — > q=oo. For the case of Eq. (|2.19p . n = —7/2. Thus, we see that the Big Rip singularity does 
not appear in finite time for the scalar field theory, although it emerges in the corresponding modified gravitational 
theory. The metric in the Einstein frame is expressed as ds 2 = e a ds 2 = —dt 2 +a (t) dx 2 , where d(t) is the scale factor 

in the scalar field theory given by a(t) = at 3 ^ n+1 ^ ^ n+2 ^ , where a is a constant. For n = —7/2, because when t — ► t s , 
t — > ^foo, the scale factor in the scalar field theory d(t) diverges at infinite time. Consequently, the (finite-time) Big 
Rip singularity in F(R) gravity is transformed into the infinite-time singularity in the scalar field theory. This shows 
the physical difference of late-time cosmological evolutions between two mathematically equivalent theories. 



III. FINITE-TIME FUTURE SINGULARITIES IN F(R) GRAVITY 

In this section, we examine several models of F(i?)-gravity with accelerating cosmological solutions ending at the 
finite-time future singularities by using the reconstruction technique explained in the preceding section. 

First, we consider the case of the Big Rip singularity [l7| . in which H behaves as H = ho/ (t s — t). Here, ho and 
t s are positive constants and H diverges at t = t s . In this case, if we neglect the contribution from the matter, the 
general solution of (|2.9[) is given by 

P(0) = P + (t s ft** + P (t s - <P) a ~ , a ± ^ - ho + 1± Wh ^±l , (3.1) 

when ho > 5 + 2^6 or ho < 5 — 2^/6 and 

P{<j>) = {t. - ^y^' 2 (At cos (( tj - 0) In ' /t " + 1 2 ° /l °' 1 ) + Bi sin ( (f> - 0) In ~^ + ^ ~ ^ ) , (3 . 2) 

when 5 + 2^6 > h > 5 - 2y/E. Here, P+, P_, A\ and B x are constants. Using Eqs. (f!Q)) , and ([2~TH)) . we find 

that when R is large, the form of F(R) is given by F{R) cx P 1 ^"-/ 2 for ho > 5 + 2^6 or /io < 5 - 2V6 case and 
F(i?) oc R( h °+ 1 ')/ 4 x (oscillating parts) for 5 + 2^6 > /i > 5 - 2^6 case. 

Next, we investigate more general singularity H ~ ho(t s —t)^ [l8| . where /io and /3 are constants, and /io is 
assumed to be positive and t < t s because it should be for the expanding universe. Even for non-integer (3 < 0, 
some derivative of H and therefore the curvature becomes singular. We should also note that in this case the scale 

factor a behaves as a ~ exp |[/io/ (ft — 1)] (ts — t)~^~^ H |, where • ■ ■ expresses the regular terms. From this 

expression, we find that if /3 could not be any integer, the value of a, and therefore the value of the metric tensor, 
would become complex number and include the imaginary part when t > t s , which is unphysical. This could tell 
that the universe could end at t = t s even if j3 could be negative or less than —1. We assume 0^1 because 
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TABLE I: Finite-time future singularities. Type I includes the case of p and p being finite at t s . In case of Type IV, higher 
derivatives of H diverges. Type IV also includes the case in which p (p) or both of them tend to some finite values while higher 
derivatives of H diverge. Here, t s is the time when a singularity appears and a 3 is the value of a(t) at t — t 3 . 



Type 




Limit 




a 


P 




IpI 




Type 


I ("Big Rip") 


t -f t s 


a 


— > oo 


P -> 


oo 


IpI- 


oo 


Type 


II ("sudden") 


t -> t s 


a 


— > a s 


P — 




IpI- 


CO 


Type 


III 


t -> t 3 


a 


— > ffls 


P — 


00 


IpI- 


00 


Type 


IV 


t -> t 3 


a 


— > ffls 


P — 





IpI- 






the case = 1 corresponds to the Big Rip singularity, which has been investigated. Furthermore, because the 
case — corresponds to dc Sitter space, which has no singularity, we take ^ 0. When > 1, the scalar 
curvature i? behaves as R ~ 12iJ 2 ~ 12/ig (f s — t)~ 213 . On the other hand, when /3 < 1, the scalar curvature i? 
behaves as i? ~ 6iJ ~ 6ho0(to — t) 13 1 . We may obtain the asymptotic solution for P when <fi — > t s : (i) For 
/3 > 1, P(<f>) ~ exp{[V2(/?-l)](i*-^ +1 }(i s -^^ 

where w = ft,o/ [2 (/3 — 1)], and A 2 and £> 2 are constants. When <f> — > i s , -P(</>) tends to vanish, (ii) For 1 > 
> 0, P(0) - B 3 cxp |- [/i /2 (1 - /?)] (t s - 4>f~ \ (i s - </-) (/3+1)/8 , where B 3 is a constant, (iii) For < 0, P(</>) ~ 

A 3 exp |- [fto/2 (1 - /?)] (t s - </') 1 " /3 } (U - ( / 1 )-( /32 - 6,3+1 ) /8 , where A 3 is a constant. Using Eqs. Q3) and (l2~TfJl) . 

we find the behavior of -F(i?) (at large i?) as summarized in Tabic II. 

In the above investigations, we found the behavior of the scalar curvature R from that of H. Conversely, we now 
consider the behavior of H from that of R. When R evolves as R ~ 6H ~ Rq (t s — t)~ 7 , if 7 > 2, which corresponds 

to = 7/2 > 1, H behaves as H ~ (t - t)" 1 ^ 2 , if 2 > 7 > 1, which corresponds tol>/? = 7-l>0,.ffis 

given by 77 ~ {Ro/ [6 (7 — 1)]} (t s — i) 7+1 , and if 7 < 1, which corresponds to — 7 — 1 < 0, we obtain H ~ if + 
{i?o/ [6 (7 — 1)]} (t s — t) 7+ , where iJo is an arbitrary constant and it does not affect the behavior of R. Hq is chosen 



to vanish in H ~ h (t s - t) p . If 7 > 2, we find a(t) oc cxp (2/7 - 1) i/Rq/U (t s - t) 11 + , when 2 > 7 > 1, a(t) 

behaves as a(t) oc exp [{Rq/ [67(7 - 1)]} (t s - i)" 7 ), and if 7 < 1, a(t) oc exp (H t + {R / [67(7- 1)]} (t s - ty 1 ^. 
In any case, there appears a sudden future singularity fis| at i = i s . 

Since the second term in if - H Q + {R a / [6 (7 - 1)]} (i s - t)~ 1+1 is smaller than the first one, we may solve Eq. (|2.9p 
asymptotically as P ~ P {l + [2h / (1 - /?)] (t s - 0) 1_/3 | with a constant P , which gives - FqR+F^ 2 ?/^^, 

where Fq and i*i are constants. When > > — 1, wc find 2/3/ (/? + 1) < 0. On the other hand, when (3 < —1, 
we obtain 2(3/ (0 + 1) > 2. As we saw in F(R) oc i? 1- "-' 2 above, for < —1, H diverges when t — > t s . Since we 
reconstruct -F(i?) so that the behavior of H could be recovered, the F(R) generates the Big Rip singularity when R 
is large. Thus, even if R is small, the F(R) generates a singularity where higher derivatives of H diverge. 

We assume that H behaves as H ~ ho (t s — t)' 13 . For > 1, when f — > t s , a ~ exp /io (t s — i) 1- ' 3 / (0 — 1) -^00 

and p ff, |p e flf| — 00. If /3 = 1, we find a ~ (t s — t) ' l0 — > 00 and p ff, |p c ff| - ► 00. If < < 1, a goes to a constant 
but p, \p\ -^00. If — 1 < < 0, a and p vanish but |p e ff| — 00. When /3 < 0, instead of H ~ /io (is — i) j on e may 
assume ~ + ho (t 8 — t) . Hence, if — 1 < < 0, p e g has a finite value 3Hq/k 2 in the limit t — * t s . If < — 1 
but /3 is not any integer, a is finite and p e s and p e ff vanish if Ho = or p e ff and p e s are finite if i?o 7^ but higher 
derivatives of H diverge. We should note that the leading behavior of the scalar curvature R does not depend on Ho 
in H ~ i/ + ho (t s — t)" 13 , and that the second term in this expression is relevant to the leading behavior of R. We 
should note, however, that Ho is relevant to the leading behavior of the effective energy density p e g and the scale 
factor a. 

In Rcf. [||, the finite-time future singularities has been classified as shown in Table I. The Type I corresponds to 
> 1 or = 1 case, Type II to — 1 < < case, Type III to < < 1 case, and Type IV to < — 1 but is not any 
integer number. Thus, we have constructed several examples of F(R) gravity showing the above finite-time future 
singularities of any type. It also follows from the reconstruction method that there appears Type I singularity for 
F(R) = R + aR n with n > 2 and Type III singularity for F(R) — R— 0R~ n with n > 0, where a and are constants. 
In fact, however, even if some specific model contains the finite-time future singularity, one can always reconstruct 
the model in the remote past in such a way that the finite-time future singularity could disappear. Positive powers of 
the curvature (polynomial structure) usually help to make the effective quintessence/phantom phase become transient 



7 



TABLE II: Summary of the behavior of F(R) gravity in case of H ~ ho (t s - t)"' 3 . Here, ci = [A /2 (/3 - 1)] (12h )~ (l3 ~ 1)n2l3) 
and c 2 = [ho/2 (1 - /3)] (-6/i /3) (/3 ~ 1)/(/3+1) . We note that -6h f3R > when ho, R > 0. 



Type I ("Big Rip") Type 


II ("sudden") 


Type III 




Type IV 


9-1 

F(7?) cx e cifl w R~i F(R) oc e 


1< p < 

-C 2 fl0+1 8(0 + 1) 


< p < 1 

(3-1 

F(i?) oce- C2R " +1 Ri 


P<- 
F(R) oc 


T, P : not integer 



and to avoid the finite-time future singularities. The corresponding examples have been examined in Refs. (3, fl8j. 

IV. CONCLUSION 

In the present article, we have reviewed finite-time future singularities in modified gravity. We have reconstructed 
an explicit model of modified gravity realizing a crossing of the phantom divide. It has been shown that the Big 
Rip singularity appears in this modified gravitational theory, whereas that the (finite-time) Big Rip singularity in the 
modified gravity is transformed into the infinite-time singularity in the corresponding scalar field theory. In addition, 
we have examined several models of modified gravity which predict accelerating cosmologies ending at the finite-time 
future singularities of all four known types. 

Acknowledgments 

The author deeply appreciates the invitation of Professor Sergei D. Odintsov to submit this article to the special 
volume Casimir effect and Cosmology on the occasion of the 70th birthday of Professor I. Brevik published by 
Tomsk State Pedagogical University, 2008. He also thanks Professor Chao-Qiang Geng, Professor Shin'ichi Nojiri 
and Professor S. D. Odintsov for their collaboration in Refs. 0, Q very much. In addition, he is grateful to Professor 
Misao Sasaki for very helpful discussion of related problems. This work is supported in part by National Tsing Hua 
University under Grant #: 97N2309F1. 



[1] S. Nojiri and S. D. Odintsov, eConf C0602061 (2006) 06 [Int. J. Geom. Meth. Mod. Phys. 4 (2007) 115] 



arXiv:hep-th/0601213 



[2] S. Nojiri and S . D. Odintsov, |arXiv:0801.4843 [astro-ph]; arXiv:0807.0685 [hep-th]; T. P. Sotiriou and V. Faraoni, 
larXiv:0805.T726l [gr-qc]; F. S. N. Lobo, larXiv:0807.T640l [gr-qc]; S. Capozziello and M. Francaviglia, Gen. Rel. Grav. 
40 (2008) 357. 

[3] M. C. B. Abdalla, S. Nojiri and S. D. Odintsov, Class. Quant. Grav. 22 (2005) L35 [arXiv:hep-th/0409177| . 
[4] L. Amendola and S. Tsujikawa, Phys. Lett. B 660 (2 008) 125. 

[5] K. Bamba, C. Q. Geng, S. Nojiri and S. D. Odintsov. la~rXiv:0810.4296l [hep-th], 

[6] S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. D 71 ( 2005) 063004 [arXiv:hep-th/0501025] . 
[7] S. Nojiri and S. D. Odintsov, Phys. Rev. D 74 (2006) 086005 [arXiv:hep-th/0608008 



[8] S. Nojiri and S. D. O dintsov, J. Phys. Conf. Ser. 66 (2007) 012005 [arXiv:hep-th/0611071| ; J. Phys. A 40 (2007) 6725 

|arXiv:hep-th/0610164| . 

[9] K. Bamba, S. Nojiri and S. D. Odintsov, JCA P 0810 (2008) 045 |arXiv:0807.2'575l [hep-th]]. 
[10] J. L. Cortes and J. Indurain. la"rXiv:0805.348"T1 [astro-ph], 

[11] S. Nojiri and S. D. Odintsov, Phys. Rev. D 77 (2008) 026007 [arXiv:0 710.1738 [hep-th]] ; G. Cognola, E. Elizalde, S. Nojiri, 
S. D. Odintsov, L. Sebastiani and S. Zerbini , ibid. 77 (2008)' 046009 arX iv:0712.40T7l [hep-th]]; G. Cognola, E. Elizalde, 
S. D. Odintsov, P. Tretyakov and S. Zerbini. larXiv:0810.4989l [gr-qc], 

[12] S. Nojiri and S. D. Odintsov, Phys. Lett. B 657 (2007) 238 [arXiv:0707.194Tl [hep-thF 



[13] F. Briscese, E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Lett. B 646 (200 7) 105 [arXiv:hep-th/0612220 
[14] S. Capozziello, S. Nojiri and S. D. Odintsov, Phys. Lett. B 634 (2006) 93 |arXiv:hep'th/0512118 



[15] Y. Fujii and K. Maeda, The Scalar- Tens or Theory of Gravitation (Cambridge University Press, Cambridge, United King- 
dom, 2003); K. i. Maeda, Phys. Rev. D 39 (1989) 3159. 



[16] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68 (2003) 123512 [arXiv:hep-th/0307288| . 

[17] B. Mclnnes, JHEP 0208 (2002) 029. 

[18] S. Nojiri and S. D. Odintsov, Phys. Rev. D 78 (2008) 046006 |arXiv:0804.35i"9l [hep-th]]. 



8 



[19] J. D. Barrow, Class. Quant. Grav. 21 (2004) L79; S. Cotsakis and I. Klaoudatou, J. Geom. Phys. 55 (2005) 306; S. Nojiri 
and S. D. Odintsov, Phys. Rev. D 70 (2004) 103522; J. D. Barrow and C. G. Tsagas, Class. Quant. Grav. 22 (2005) 1563; 
M. P. Dabrowski, Phys. Rev. D 71 (2005) 103505; Phys. Lett. B 625 (2005) 184; L. Fernandez- Jambrina and R. Lazkoz, 
Phys. Rev. D 70 (2004) 121503; 74 (2006) 064030; P. Tretyakov, A. Toporensky, Y. Shtanov and V. Sahni, Class. Quant. 
Grav. 23 (2006) 3259; M. Bouhmadi-Lopez, P. F. Gonzalez-Diaz and P. Martin-Moruno, Phys. Lett. B 659 (2008) 1; 
C. Cattoen and M. Visser, Class. Quant. Grav. 22 (2005) 4913; I. Brevik and O. Gorbunova, Eur. Phys. J. C 56 (2008) 
425; J. D. Barrow, Phys. Lett. B 235 (1990) 40. 



